First Principle Noncollinear Transport Calculation and Interfacial Spin-flipping of 

Cu/Co Multilayers 
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In this paper the first principle noncollinear transport calculation for Cu/Co(lll) including in- 
terfacial spin-flipping was performed. We modeled spin-flipping at the interface by assuming a 
noncollinear magnetic structure with random magnetization orientation which satisfied Gaussian 
distribution along average magnetization direction. The relationship between spin-dependent con- 
ductance including interfacial spin-flipping and random magnetization orientation distribution width 
was obtained. For certain distribution width, our defined spin-flipping ratio coincides with the range 
of experimental spin-flipping probability P = 1 — e~ 5 , where 8 = 0.25 ± 0.1. The magnetoresistance 
in Co/Cu/Co spin valve system including interfacial spin-flipping has also been calculated. 
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I. INTRODUCTION 

The electron transport across ferromag- 
netic/nonmagnetic (FM/NM) interface is of inter- 
esting in the past two decadesi^ Based on the two 
current model(non spin-flipping theory), ab initio 
calculation with no-free parameter of interfacial specific 
resistance^ agreed reasonably well with experimental 
data for some lattice-matched metal pairs, such as fee 
Cu(lll)/Co(lll)£ and bec Fe(110)/Cr(110)£ 

The spin-flipping at the FM/NM interface, which 
certainly exists in the experiment, is not well stud- 
ied, partially due to the theoretical difficult, and par- 
tially due to lacking of reliable experimental data J How- 
ever, spin-flipping at the FM/NM interface is getting 
increasing importance. Spin-flipping at the ferromag- 
netic/supperconductor interface can induce spin triplet 
pairing in the ferromagnetic sidei2i2ii£ In addition, the 
spin- flipping at the FM/NM interface can also change 
the spin torques induced by current ] 11 ^ 2 ' 13 

Moreover, Geux et al. have calculated the trans- 
mission probability in the presence of magnetic impu- 
rity scattering with spin-flipping by the effective mass 
approximation.— They calculated the transmission prob- 
ability matrix in spin space and obtained that the spin- 
flipping probability is proportional to the impurity den- 
sity. They also found that to the first order the calcu- 
lated conductances from transmission probability matrix 
decrease linearly with increasing the impurity density. 
However, the transport properties through the real inter- 
face with interfacial spin-flipping process have not been 
studied yet. 

In this paper, we will calculate the scattering matrix 
of real FM/NM interface with spin-flipping process by 
the first principle noncollinear transport calculation. 15 
The spin-flipping at the interface is modeled by assum- 
ing a noncollinear magnetic structure with random mag- 
netization orientation which satisfied Gaussian distri- 
bution along average magnetization direction. By the 
noncollinear transport calculation,^ we can obtain how 
the interfacial conductance changes with random interfa- 
cial magnetic structure and the effect of interfacial spin- 



flipping on the magnetoresistance in FM/NM/FM spin 
valve system. 



II. COMPUTATIONAL DETAILS 

Our calculation of scattering matrix is based on the 
surface Green's function method 16 with tight-binding lin- 
ear muffin tin orbital basisJ^ First, the self-consistent 
one-electron effective potential in our calculation is ob- 
tained from collinear electron structure calculation 16 
without spin-orbit coupling. The Hamiltonian Hq is con- 
structed by this self-consistent potential and is diagonal 
in spin space. Second, the rigid potential approxima- 
tion has been used in our noncollinear transport calcula- 
tion. In this approximation, we rotate the Hamiltonian 
Hq which is in local quantum axis representation in spin 
space to the global quantum axis representation. So the 
Hamiltonian for transport calculation can be written as 
H' = U(9,ip)H Ui(9,<p), where U(6,<p) is the unity ro- 
tation matrix in spin space. 6 and ip are the polar angle 
and azimuth angle of the local quantum axis respectively 
(global quantum axis is taken as z axis and tp = in 
our calculations). Therefore the spin-flipping is only in- 
duced by interfacial magnetic disorder and the spin-orbit 
coupling is neglected in our transport calculation. 

Due to the lack of the details of interfacial magnetic 
disorder, we assume that spin-flipping is introduced by 
the Gaussian random distribution of magnetization ori- 
entation. In our calculation all the magnetization orien- 
tation of Co atoms are along one global quantum axis ex- 
cept the Co monolayer at the interface. Further the disor- 
dered magnetization orientation is modeled by 10 x 10 lat- 
eral supercell. The deviation angles within the 10 x 10 su- 
percell satisfy the Gaussian random distribution, where 
the average orientation is along global quantum axis and 
the distribution width is AO. So the most deviation an- 
gles in the supercell are in the range of (-A6, Ad). How- 
ever, the magnitude of magnetization for each Co atom 
in the supercell is constant. In this paper we have calcu- 
lated the transport properties for different magnetization 
orientation distribution width, which is from AO = to 
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(a) sharp interface 
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(c) roughness interface^ 
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(b) sharp interface 
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(d) roughness interface 
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FIG. 1: The spin-dependent conductances of sharp and roughness Cu/Co(lll) interfaces, (a) and (c) The conductances will 
be saturated for large A9 case and the critical A9 of roughness interface is larger than that of sharp interface, (b) and (d) For 
both sharp and roughness interface the spin in majority channel can be flipped more easily than that in minority channel. 



180 (deg). Here the larger distribution width A8 implies 
the larger probability of spin-flipping scattering process 
at interface. 

The noncollinear magnetic structure at interface has 
been calculated by other group using first principle lo- 
cal spin density calculations,— and the multiple and 
metastable noncollinear magnetic structure have been 
obtained in copper-permalloy interface. There are four 
metastable noncollinear states in which the energies are 
lower than the energy of collinear state or total random 
state. The collinearity of those four metastable states in 
their paper— is in the range of about 0.4 to 0.8, which 
corresponds to our parameter A8 as 40 to 70 (deg). 



III. RESULTS AND DISCUSSION 

In our noncollinear Cu/Co(lll) transport calculation 
the lattice constant is taken as a=3. 549A. The spin- 



dependent conductance is 

where i a ° is transmission matrix element for bloch 
state (v, a') in lead Cu to bloch state (/i, a) in lead 
Co and k|| is lateral wave vector. FigureQ] shows the 
spin-dependent conductances with different magnetiza- 
tion orientation distribution width A9, where the rough- 
ness interface is modeled by 2ML of 50%-50% alloy 
in a 10x10 lateral supercell, which can be denoted as 
Cu[Cuo.5Coo.5|Cuo.5Coo.5]Co. For this roughness inter- 
face the random Gaussian distribution of magnetization 
direction only takes place at the interfacial magnetic 
atoms (Co atoms). Here G'" and are conductances 
for majority and minority electron channel with unchang- 
ing the orientation of spin. G^(G^) is conductance for 
spin-flipping process, which describes the probability of 
majority (minority) being scattered to minority (major- 

ity). 
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FIG. 2: The spin-flipping ratio across Cu/Co(lll). The solid 
circle is the spin- flipping probability value P = 1 — e~ s m 1 — 
g-0-25 ^ g 22 inferred from experiment data and the error bar 
comes from the uncertainty S — 0.25±0.1. The solid star is the 
estimated spin-flipping probability only due to the spin-orbit 
coupling. The range of AO between dash line corresponds 
to the collinearity 0.4 to 0.8 which is result of the magnetic 
structure calculation. 



As shown in FigureQJa) and (c), for both sharp and 
roughness interfaces, with increasing the magnetization 
direction distribution width AO the spin-flipping conduc- 
tances increase and the majority (minority) conductance 
decreases, where the majority conductance decreases 
more rapidly and the conductance increases also 
more rapidly than Gr^. For sharp interface with increas- 
ing AO, the total conductance (G iT + G TT + G n + G u ) 
increases slightly at first and decreases to constant when 
AO > 90 (deg). But for the roughness interface, the 
total conductance decreases monotonically and the crit- 
ical distribution width AO where the conductance start 
to be saturated is about 120 (deg), which is larger than 
that for sharp interface. In our results of electron struc- 
ture calculation, the average interfacial magnetic moment 
per atomic sphere is about 0.77/xs for roughness inter- 
face and 1.58 \ib for sharp interface. Moreover for the 
same distribution of random interfacial magnetic struc- 
ture, the larger interfacial magnetic moment will lead 
to more strongly spin-flipping scattering of the incoming 
electron. Therefore, it needs much more degree of inter- 
facial magnetic structure disorder to saturate the con- 
ductance in the roughness interface. 

In the case of non spin-flipping Cu/Co(lll) interface 
the minority electron is reflected more strongly than the 
electron in majority channel. Taking into account the 
scattering induced by noncollinear magnetization, the 
conductances G^* and G^ both decrease with increas- 
ing AO. Further the influence of noncollinear magnetiza- 
tion scattering is relatively more important for majority 
channel than for minority channel. Therefore as shown 



FIG. 3: The rounded interfacial specific resistance 2ARp/M 
for different distribution width AO, where the line in the figure 
corresponds to the value from the experiment. The region 
between two dash line corresponds to the result of magnetic 
structure calculation. 



in FigureQJa) and (c), G'* decreases more rapidly than 
G^. In addition, for both sharp and roughness interfaces 
we observe that G iT /G TT > G u /G u in FigureQJb) and 
(d), which also indicates that the spin in majority channel 
can be flipped more easily than that in minority channel. 

The interfacial spin-flipping parameter r5 usually de- 
scribes the interfacial spin-memory-loss which is defined 
as 5 = ti/l s f within VF theory— and the spin-flipping 

probability is^ P — 1 — e~ s . For Cu/Co interface 5 ~ 0.25 
which is indirectly inferred by explaining the difference 
between 'interleaved' and 'separated' sample within VF 
theoryiSSiSLSS. There is no direct measurement of pa- 
rameter S for F/N interface up to date and the value 
of 5 inferred from experiment is quite uncertain^ In 
this paper, we defined the spin-flipping ratio as P = 
(G iT + G Ti )/(G TT + G u ). From Figured one can see 
that for AO « 40 to 70 (deg) our defined spin-flipping 
ratio coincides with the range of spin-flipping proba- 
bility inferred from experimental data (P = 1 — e~ s , 
S = 0.25 ± 0.1). 20 This might be just an accident coin- 
cidence because we have not calculated the spin-flipping 
probability directly and our definition of the spin flipping 
ratio is only used for giving some quantitative informa- 
tion of spin-flipping. In addition, K. Eid et a/.— esti- 
mated the interfacial spin-flipping probability at Cu/Co 
interface only due to the spin-orbit coupling and the re- 
sult is P = 1 - e~ s w 1 - e-°- 2 w 0.18, which is on the 
small side of the experimental P value range. The esti- 
mated spin-flipping probability due to spin-orbit coupling 
approximates the P value for the case of AO = 40 (deg) 
in our calculation and is much smaller than the satu- 
rated spin-flipping probability due to magnetic disorder. 
Considering our calculation only taking the interfacial 
magnetic disorder into account, it is suggested that the 
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FIG. 4: The calculated MR of the spin valve system 
Co/Cu/Co with different interfacial spin-flipping, where the 
two Cu/Cu interfaces are both roughness with 2ML 50%-50% 
alloy. The region between two dash line also corresponds to 
the result of magnetic structure calculation. 



spin mixing effect of interfacial magnetic disorder is more 
prominent than that of spin-orbit coupling on the Cu/Co 
transport properties, e.g. current-perpcndicular-to-planc 
(CPP) magnetoresistance (MR). 

Figure[3] shows the rounded interfacial specific resis- 
tance which is defined as^ 
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where A is section area, Nq u , Ng a and Nq are 
the Sharvin conductances. First, in non spin-flipping 
case for sharp interface we obtain 2ARp/^ ~ 0.96 
(lCP 15 f2TO 2 ) which is close to the experimental valued ~ 
1.0 (10 _15 f2m 2 ). Second, due to the additional scattering 
from noncollinear magnetization, with increasing A9 the 
rounded interfacial specific resistance increases and for 
large distribution width the calculated 2AR F / N w 1.23 
(10~ 15 f2m 2 ) which is about 128% of non spin-flipping 
specific resistance. Moreover the experimental specific 
resistance value corresponds to the case of A9 « 40 (deg) 
in our calculation, which indicates that the spin-flipping 
at interface can also explain the deviation of rounded 
interfacial specific resistance between sharp interface cal- 



culation and experimental data. The rounded interfacial 
specific resistance of roughness interface is also shown in 
FigurelU It can be seen that the specific resistance of 
roughness interface also increases with increasing A9 as 
the case of sharp interface. However, the value of rough- 
ness interface is larger than that of sharp interface and 
the A9 in which the specific resistance start to be satu- 
rated is also larger than that of sharp interface. 

Figured] shows the magnetoresistance(MR) of the spin 
valve system Co/Cu/Co with roughness interface for 
different interfacial spin-flipping. We have calculated 
the conductances of the parallel G(P) and antiparallel 
G(AP) configuration, and the MR is defined as MR = 

G(P G(AP) AP) x 100% - 0ne can observe that the MR de- 
creases rapidly with increasing distribution width A9 and 
for A9 > 90 (deg) the MR is nearly constant about 5%. 
For the region of A9 — 40 to 70 (deg) corresponding to 
the result of magnetic structure calculation^ the MR is 
in range of 11.7% to 6.7%, which is about 2/3 to 1/3 of 
the non spin-flipping MR ~ 18.8%. 



IV. CONCLUSION 

In this paper the first principle noncollinear transport 
calculation for Cu/Co(lll) including interfacial spin- 
flipping was performed. We modeled spin-flipping at the 
interface by assuming a noncollinear magnetic structure 
with random magnetization orientation which satisfied 
Gaussian distribution along average magnetization di- 
rection. The relationship between spin-dependent con- 
ductance including spin-flipping conductance and ran- 
dom magnetization orientation distribution width was 
obtained. We found that the conductances start to be 
saturated when the distribution width A9 larger than 
the critical width. For distribution width Af5=40 to 70 
(deg) our defined spin-flipping ratio coincides with the 



range of spin-flipping probability P 



1 



where 



5 = 0.25 ± 0.1 inferred from experimental data. In 
addition we also found that the magnetoresistance of 
Co/Cu/Co spin valve system decreases rapidly with in- 
creasing interfacial spin-flipping probability. 
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